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ABSTRACT 

We report high-cadence time-series photometry of the recently-discovered transiting exoplanet sys- 
tem HD 17156, spanning the time of transit on UT 2007 October 1, from three separate obser- 
vatories. We present a joint analysis of our photometry, previously published radial velocity mea- 
surements, and times of transit center for 3 additional events. Adopting the spectroscopically- 
determined values and uncertainties for the stellar mass and radius, we estimate a planet radius 
of Rp = 1.01 ± 0.09i?jup and an inclination of i = 86.5^q y degrees. We find a time of transit center 
of Tc = 2454374.8338 ± 0.0020 HJD and an orbital period of P = 21.21691 ± 0.00071 days, and 
note that the 4 transits reported to date show no sign of timing variations that would indicate the 
presence of a third body in the system. Our results do not preclude the existence of a secondary 
eclipse, but imply there is only a 9.2% chance for this to be present, and an even lower probability 
(6.9%) that the secondary eclipse would be a non-grazing event. Due to its eccentric orbit and long 
period, HD 17156b is a fascinating object for the study of the dynamics of exoplanet atmospheres. 
To aid such future studies, we present theoretical light curves for the variable infrared emission from 
the visible hemisphere of the planet throughout its orbit. 

Subject headings: planetary systems - stars: individual (HD 17156) - techniques: photometric 



1. INTRODUCTION 

Twenty-eight transiting planets are now reported in 
the literature^, and the doubling time scale for new 
detection s is now roughly one year. As reviewed by 
[Charbon neau et al.l (j2007l ). it is these objects that have 
allowed the study of the physical structure and atmo- 
spheric chemistry and dynamics of gas giant exoplan- 
ets, and opened the field of comparative exoplanetology. 
However, the majority have orbital periods of a few days, 
due both to the lower geometric probability of transits 
occurring as the orbital period increases, and due to the 
limitations of ground-based photometric transit surveys 
that strongly favor detection of short-period systems. 
The latter problem can be circumvented by searching 
for transits of known radial velocity planet-bearing stars, 
putting to use the radial velocity information to constrain 
the possible time of transit. Th is is the primary aim of 
the Transitsearch.org network (jSeagroves et al.i [20031 )^. 



which, starting in 2002, has been conducting photomet- 
ric follow-up observations of known radial velocity de- 
tected planets. This network is a prime example of suc- 
cessful collaboration between amateur and professional 
astronomers. 

Although longer period planets possess a lower geomet- 
ric probability to transit, there exists a loop-hole for some 
planets on highly eccentric orbits. The transit probabil- 
ity for a planet on an eccentric orbit, with periastron 
near inferior conjunction, is amplified by a factor 



A = 



1 + ecos(^ — oj) 



(1) 



Electronic address: |jirwin at of a dot harvard dot edu] 
^ Alfred P. Sloan Research Fellow 
^ See http : //www. inscience . c h/translts/ | 
^ See http : //www. t ransitsearch. org| 



where e is the eccentricity, and uj is the argument of peri- 
center (Seagroves et al. 2003). An amplification factor 
near 2.9 for the planet orbiting HD 17156 (Fischer et al.l 
|2007( ). along with an extremely short time window for 
possible transits, brought the system to the attention 
of the Transitsearch.org network. This was rewarded 
by the detecti on of transits as recently announced by 
iBarbieri et all (|2607D . 

This unique system contains a 3.1-Mjup transiting 
planet in a 21.2 day highly-eccentric orbit (e = 0.67) 
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about a bright {V = 8.2) GOV star. iGillon et al.l (pOOTh 
have recently presented refined estimates of the system 
parameters based on new photometric data. Even among 
planets with periods beyond 5 days, where the timescale 
for tidal circularization quickly grows to exceed the age 
of most systems, HD 17156b's eccentricity is unusually 
high, thus making it an interesting test case for models of 
planetary migration. Preliminary measurements of the 
Rossiter- McLaughlin effect (Rossitcr 1924; McLaughlin 
fl92l iGaudi fc WinnI [2007L fWinn 2003) bv lNarita et all 
(I2007D show evidence for a misalignment between the 
stellar spin axis and the planetary orbit axis, which may 
indicate a migratio n mechanism involving planet-planet 
scattering (see e.g. IChatteriee et al.|[2007[ ). or Kozai mi- 
gration under perturbation by a yet undetected stel- 
lar companion or second planet (iFabrvckv fc Tremain"i 
l2007HWu. Murrav fc Ramsahaill2007D . The large eccen- 
tricity and long period also make HD 17156 a particularly 
attractive target for several additional follow-up studies. 
First, it presents a unique opportunity for the study of 
the structure and dynamics of a gas-giant atmosphere un- 
der strongly varying illumination, thro ugh infrared mon- 
itorin g of the planetary emission (e.g. | Harrington et al.l 
[200l lCowan et al.l[2007t iKnutson et al.l l2007al lbf) . Sec- 
ond, by monitoring successive times of transit, the pres- 
ence of additional bodies in the system can be detected or 
constrained from t he resulting perturbation s on the or- 
bit o f HP 17156b (iHolman fc Murray! 120051: lAgol et all 
\2(M iSteffen fc Ago]ll2005l i! 

The purpose of our paper is three-fold: First, we 
seek to refine the estimates of the system parameters 
that were only poorly constrained by the light curves of 
iBarbieri et all (|2007f ). Second, we document the time 
of center of transit and search for transit timing varia- 
tions. Third, we evaluate the likelihood that a secondary 
eclipse will be observable for HD 17156, and present the- 
oretical predictions of the infrared phase variations as 
might be detected with the Spitzer Space Telescope. We 
elected to perform an independent analysis from that of 
iGillon et al.l (|2007f l. and therefore do not use their orbital 
parameters as constraints, since at the time of writing 
these are still subject to change as their paper has not 
yet been accepted for publication. We do however make 
use of the time of mid-transit presented in their Table 1. 

We begin by presenting a description of our observa- 
tions of the transit event on UT 2007 October 1 in f|2] We 
then present our combined analysis of the extant radial 
velocity measurements (Sj3]) and photometric data (Sj4]), 
and our resulting estimates of the system parameters and 
likelihood of a secondary eclipse. In ^J5l we compare our 
constraints on the planet mass and radius with planetary 
structural models, summarize the current constraints on 
the presence of transit timing variations in the system, 
and present model predictions for the planetary infrared 
light curve. We conclude in ^with a discussion of com- 
pelling avenues for future research. 

2. OBSERVATIONS 

Using the orbital period given bv lFischer et al.l (120071) 
and th e time of mid-transit measured bv IBarbieri et all 
()2007t ). we predicted that an event would occur on UT 
2007 October 1 with transit center at UT 7:53, which 
was well-situated for observatories in the southwestern 
United States. We present below a description of the 



data we gathered from three such observatories, and sum- 
marize the calibration of the raw frames and extraction 
of the time series for each. 

2.1. Mount Laguna Observatory 

We observed HD 17156 with the Smith 0.6 m tele- 
scope at Mount Laguna Observatory using a Bessell R 
filter and an SBIG STL-IOOIE CCD camera (giving a 
scale of O'.'4/pixel at the //20 Cassegrain focus), using 
BD+71 168 {V = 9.57, B8 spectral type) as our local 
comparison star. The integration time for each exposure 
was 4 s, with a 4.3s readout time and a net cadence of 
8.3 s. For each frame, we calculated the time at mid- 
exposure and corrected this time to heliocentric Julian 
Day. We dark-subtracted and flat-fielded the raw images 
in IRAF, and then performed aperture photometry using 
a circular aperture with radius 12.8 arcsec. The use of 
the large aperture insulated our analysis against seeing- 
induced aperture losses, as the seeing varied from 2.1 arc- 
sec to 2.7 arcsec over the course of the night. Clouds 
were present at several times during the observations, 
especially after mid-transit. To avoid systematic effects 
due to cloud-induced differential extinction of the differ- 
ent color stars, we rejected observations for which the 
calibration star showed more than 8% light loss. The 
resulting 2148 points were binned 20:1 and uncertainties 
were estimated from the RMS scatter in each bin divided 
by '\/20. Typical uncertainties for the binned points were 
0.3-0.7%. 

2.2. Torrance California 

We observed a field centered on HD 17156 using a com- 
mercial 0.25m Meade LX200 telescope and a Meade Deep 
Sky Imager Pro II CCD camera. We employed an f/3.3 
focal reducer to increase the field-of-view to ~ 20 x 30 ar- 
cmin. We gathered 10 345 exposures through an Edmund 
Optics IR-pass filter (approximating the conventional I- 
band), between UT 4:20 and 13:25. The integration time 
for each image was 2 s. Weather conditions were clear. 
Approximately 10% of the frames were lost due to an in- 
termittent issue with the CCD readout, and due to wind 
shake, and the guide star was lost temporarily at 6:10 
UT, resulting in a large pointing shift (~ 5 arcmin). 

These i mages were proce ssed using the software de- 
scribed inl lrwin et al.l (j2007f ). which was or iginally devel- 
oped for the Monitor project (jAigrain et al. 2007). In- 
dividual frames were processed by subtracting a master 
dark and dividing by a master flat fleld . We then ran the 
pipeline source detection software (see llrwin "eTalllMl 
on all of the data frames, choosing a single frame taken 
in good conditions for use as a reference. For each image, 
we derived an astrometric solution employing an imple- 
menta tion of the triangle matching algorithm of iGrothI 
()1986[ ) to cope with shifts in the telescope position (these 
were ~ 10 pixels peak-to-peak, in a periodic fashion, pre- 
sumably due to errors in the telescope worm drives) , and 
fleld rotation due to polar misalignment of the mount. 

We generated light curves using the method described 
in llrwin et al.l (j2007f ). using the nearby star BD-l-71 168 
as the comparison source for the differential photometry. 
We found that the normalized time series data gathered 
before and after the large image motion at UT 6:10 have 
an photometric offset of roughly 3%. Since the data ob- 
tained before UT 6:10 were well before the ingress, we 
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simply trimmed them from the time-series; similarly, we 
trimmed the data that occurred more than 0.2 days af- 
ter the transit midpoint. We also trimmed single large 
outliers by identifying all points that deviated from the 
median of the time-series by more than 8%. The typi- 
cal photometric errors in the final, trimmed time-series 
are 2% per data point, as determined from the RMS 
variation from the data gathered after egress. We then 
binned this time series by a factor of 15 (and assumed 
that the errors in the binned data were reduced accord- 
ing to Poisson statistics); this binning was necessary to 
keep the computational time for the model fitting (§4|) to 
a manageable duration, but the binning did not degrade 
the quality of the fits. 

2.3. Las Cumbres Observatory Global Telescope 

We also observed the event using a 0.4m Meade 
RCX400 telescope on a custom mount, temporarily lo- 
cated in the car park of the Las Cumbres Observatory 
Global Telescope (LCOGT) offices in Goleta, California. 
We used an SBIG STL-6303E CCD camera to image a 
20 X 30 arcmin field with 0.6-arcsec pixels. 380 useful 
images were taken between UT 5:30 and 11:00, with the 
airmass decreasing from 1.6 to a minimum of 1.25 near 
UT 10:00. No autoguiding was used, and the field drifted 
by ^ 100 arcsec during the run. The observations were 
made through the "Red" filter from the LRGBC filter 
set provided with the camera by SBIG,^ with ^ 90% 
transmission in the range 580-680 nm. Exposure times 
varied from 20 to 30 seconds. The readout time was 22 
seconds. The telescope was defocused to avoid saturat- 
ing the target star, which was the brightest in the field. 
The amount of defocus was increased at about UT 06:30, 
shortly before the transit ingress, but this did not signif- 
icantly affect the relative photometry at that point. 

Bias and dark subtraction and fiat fielding were done 
using standard IRAF tasks. HD 17156 and ten refer- 
ence stars were measured using the DAOPHOT aperture 
photometry package, with apertures of radius 10.6 arc- 
sec. The light curve was divided by the average of the 
reference light curves to remove the effects of varying at- 
mospheric extinction and focus changes. The standard 
deviation of the (unbinned) points outside the transit was 
0.4%. 

The scatter of these data is lower than for the Mount 
Laguna 0.6m observations, despite the smaller aperture, 
since the larger field-of-view permits the use of multi- 
ple comparison stars. In the Mount Laguna data, the 
signal-to-noise ratio of the final light curve was limited 
by that of the comparison light curve, due to the use of 
a comparison star fainter than the target. By using mul- 
tiple comparison stars for the Las Cumbres data, it was 
possible to attain a higher signal level in the comparison 
light curve, such that the overall noise was dominated 
by the target, hence giving a net overall improvement 
in signal-to-noise ratio despite the use of a smaller aper- 
ture telescope. Observing conditions were also somewhat 
worse at Mount Laguna, and this along with the use of 
defocus to render the observations insensitive to seeing 
variations, and to average out pixel-to-pixel flat fielding 
errors, also contribute to reducing the noise level in the 
Las Cumbres data. 

* |http: //www. sblg. com/large_f ormat/f llterchart_large .htm| 



3. ANALYSIS OF THE RADIAL VELOCITY DATA 

Our orbital model was obtained by jointly fitting the 
iFischer et al.l ([2007) radial velocity data (converted to 
HJD) in conjunction with the following four observed 
central transit times, numbered N = —1,0,2,3 (where 
we denote our observation on UT 2007 Octobe r 1 as event 
N = 0): Tc_i = HJD 2454353.61 ±0.02 (Bar bieri et all 
l2007l) . Tc.o = HJD 2454374.8338 ± 0.0020 ( this paper; 
see SI), Tc,2 = HJD 2454417.2645 ±0.0021 (,Narita et all 
[20071) , andTc.3 = HJD 2454438.4835l[! [J[5^^ (|Gillon et all 
[2001 . 

For a Keplerian orbit, the instantaneous stellar ra- 
dial velocity is given by Vmod.i — K [cos(/i + w) + 
e cosw], where fi is the true anomaly of the planet 
at time ti. If we assume an edge-on orbit, the time 
of central transit occurs when ft = 7r/2 — w. Given 
the precision of the available photometry, this approx- 
imation to ft is excellent. The true anomaly, /, is 
in turn related to the the eccentric anomaly, i?, via 
E = 2 arctan v'(l - e)/(l + e) tan(//2), so that a cen- 
tral transit, Tc.i^ occurs at a fixed interval AT = P{E — 
esin£^)/27r following the epoch of a periastron passage. 



p,2 ■ 



Given Njiv radial velocity measurements (here, 
N]:(y — 33) and Nt^ central transit measurements, the 
goodness-of-fit function of a particular orbital model is 
calculated as: 



Nrv /-t^ 

( ^mod,j: 

1=1 



-E 

1=0 



T 



AT -Tr 



We use a Levenberg-Marquardt algorithm to minimize 
X^. The best-fitting orbital parameters are listed in Ta- 
ble [TJ To obtain the quoted uncertainties, a simple boot- 
strap procedure was used. An aggregate of 1000 alter- 
nate datasets was created by (1) redrawing the radial 
velocities with replacement, and (2) sampling T^ values 
for the four transits by drawing from the Gaussian dis- 
tributions implied by the quoted errors. Fits to each 
of these data sets were obtained, and the resulting dis- 
tributions of orbital parameters were used to derive the 
68.3% confidence intervals for each parameter. In Table 
(U we have also given in brackets the standard deviations 
of these parameter distributions. The latter are larger 
than the former due to the existence of outliers in the 
bootstrap sample. These are generated because the rel- 
ative importance of each radial velocity data point in 
constraining the derived orbital parameters depends on 
the orbital phase as a result of the highly eccentric orbit: 
the points close to periastron have the largest effect on 
the fit. Consequently, when the bootstrapping procedure 
exchanges a point close to periastron, the parameters es- 
timated from this particular realization of the procedure 
show sizable variations. We therefore prefer the 68.3% 
confidence intervals, since these are robust to the pres- 
ence of the outliers, which arc in turn the result of the 
bootstrap procedure not being strictly applicable to data 
such as these which are not independent and identically- 
distributed. 

We note that the addition of the transit timing data 
provides a very strong constraint on the orbital period, 
but has only a minor effect on the other orbital parame- 
ters, which are in excellent agreement with the solution 
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TABLE 1 
Fitted Parameters for HD 17156. 



Parameter^ Value 



P 21.21691 ± 0.00071 (0.00094) d 

e 0.670 ± 0.006 (0.020) 

£j 121?3±0?9 (1?8) 

K 273.8 ± 3.4 (16.3) ms'i 

Tj, 2453738.605 ± 0.024 (0.036) HJD 

Tc 2454374.8338 ± 0.0020 HJD 

i 86?5+J-^^ 

Rp/Ri, 0.070 ± 0'.003 

Rp 1.01 ± 0.09 Rjup 

Ri, 1.47 ±0.08^0 

P 3.8l?;f g cm-3 

/i 0.9989 ± 0.0025 

fei -0.0006 ± 0.0015 

/2 1.006 ±0.003 

fea 0.0049 ± 0.0022 

Cl -0.0024 ± 0.0005 

C2 0.0173 ± 0.0029 

C3 0.204 ±0.064 



For the quantities P, e, w, and Tp, 
68.3% confidence intervals are quoted, with 
the standard deviation from bootstrapping 
in brackets. The latter are larger due to 
outliers in the bootstrap sample generated 
by the high sensitivity of the fitted param- 
eters to the velocities near periastron pas- 
sage. As discussed in the text, the 68.3% 
confidence intervals are preferred. 

given bv lFischer eTal] (|2007f) . 

4. PHOTOMETRIC ANALYSIS 

We pooled together the three photometric data sets to 
fit a model light curve parametrized by the orbital and 
physical properties of the HD 17156 system. Our model 
incorporates 4 orbital parameters initially determined 
through the fit to the RV data: the period P, eccentric- 
ity e, time of periastron Tp, and argument of pericenter 
Lo. Given this orbit, we calculate a baseline light curve 
for a quadratically limb-darkened star (characterized by 
the stellar mass M^,, radius i?*, and 2 limb darkening co- 
efficients for each photometric band pass) that is being 
occulted by a planet (rp = Rp/Ri,) orbiting at an incli- 
nation i. To calculate t he baseline light curve, we employ 
the analytic formulae of lMandel fc Agoll (|2002D , together 
wi th the quadratic lirn b-darkening coefficients tabulated 
bv lClaretl ((2000l[200l . for Toff = 6000 K, [Fe/H] = 0.2, 
logfl = 4.0 (based on the stellar parameters given by 
iTakeda et al.ll2007l ). We adopted i?-band coefficients for 
the Mount Laguna data, SDSS r-band coefficients for 
the Las Cumbres data, and J-band coefficients for the 
Torrance California data. In the latter two cases, the 
filters used for the observations only approximate r and 
/, but the approximation is more than suitable for our 
purposes (in particular, the goodness-of-fit statistic de- 
scribed below is negligibly affected by errors in the as- 
sumed bandpass). The baseline flux is then modified by 
observational correction factors unique to each dataset 
(7 additional parameters described below). 

Nominally, our model fixes P, e, lo, 6 limb darkening 
coefficients, and M^,, although we iteratively update the 



first three of these parameters as the RV fit is apprised 
of the transit timing resulting from the light curve fit. 
Ultimately, the change in these three parameters from 
the iterative update process had negligible effect on the 
quality of the light curve fit and had negligible conse- 
quences for the stellar and planetary properties deter- 
mined by the transit analysis. Note that although Tp 
is a fit- variable in our transit analysis, the value that is 
really being constrained by the photometry is the time 
of central transit, which is related to a degenerate 
combin ation of P, e, cu, and Tp. We fix to the value 
1.2Mq (jFischer et al.ll2007l ). which was determined by 
matchi ng spectroscopic ob servations to stellar evolution 
tracks (iTakeda et al.ll2007l) . Here, the uncertainty in M^,, 
for which the 95% confidence interval is 1.1 — 1.3 Mq, has 
only weak impact on the light curve fit. This is because, 
as described below, we adopt an external constraint on 
the stellar radius. 

Our light curve model employed 11 free parameters: 
Tp, the planet-star radius ratio rp, P^,, i, and 7 addi- 
tional parameters related to observational corrections. 
The Torrance and Mount Laguna fluxes were each mul- 
tiplied by airmass corrections of the form / x exp(— fca), 
where / is a normalization factor, a is the airmass, and 
k is the extinction coefficient. We note, however, that 
including the Mount Laguna airmass correction had lit- 
tle effect on the quality of fit (indeed, the determined 
/ and k were consistent with 1 and respectively; see 
Table [T]). We attempted the same airmass correction 
for the Las Cumbres data, but found that a quadratic 
function of time significantly improved the quality of fit. 
We adopted a function as our goodness-of-fit statistic, 
with an additional term reflecting a spectroscopic prior 
on P*, which is approximately Gaussian with tr — 0.085 
(G. Takeda, personal communication). This additional 
constraint proved necessary, as our data could not inde- 
pendently determine the stellar radius. The goodness-of- 
fit function is 

2_^f fmod{t)~fobs{^) y (P./P0 - 1.47)^ 

^ ^{ a, J ^ 0.0852 ^-^l 

where jmodii) is the calculated flux at the time of the i*'* 
data point, /ofcs(*) is the i^^ flux measurement, and oi is 
the uncertainty of the i*'' flux measurement (derived as 
described for each data-set in [J2])- Using an IDL imple- 
ment ation of the cunoeba algorithm (e.g. see lPress et alJ 
Il992t l. we performed an initial minimization over the 
space of free parameters in our model. Using the results, 
we then rescaled the Oi so that the reduced equals 1 
separately for each dataset. We performed iterations of 
transit model fitting in conjunction with radial velocity 
fitting, using the derived transit timing to inform the ra- 
dial velocity fit, which in turn updated the P, e, and w 
used for the transit fit. As mentioned above, the param- 
eter update had negligible effect on the derived stellar 
and planetary properties. The data, corrected for the 
airmass and instrumental effects, are shown in Figure [T] 
along with the best-fitting solution. 

The uncertainty in the transit time was assessed by 
finely stepping Tp through values near its best fitting 
value and calculating the minimum each step, al- 

lowing the other free parameters to vary. The 1-cr upper 
and lower limits were assessed by noting the Tp values at 
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Los Cumbres 0.4m 



Torrance CA 0.25rT 



-0.2 -0.1 0.0 0.1 0.2 

time from center of tronsit (days) 

Fig. 1. — Binned light curve of HD 17156, with the airmass 
and instrumental corrections described in §4 applied. Top: Mount 
Laguna ij-band photometry. Middle: Las Cumbres photometry in 
a filter approximating r. Bottom: Torrance CA photometry in a 
filter approximating I. Dashed lines show the model fit simulta- 
neously to the light curve and radial velocity data, as discussed in 
the text. 

which Ax^ = 1- Given the other orbital parameters, we 
mapped each Tp into the corresponding time of central 
transit, Tc- The resulting best fitting Tc and uncertain- 
ties are reported in Table [TJ Over the range of allowed 
Tp, we found very little correlation between it and the 
best fitting stellar/planetary parameters. For this rea- 
son, we decided to fix Tp at its best fitting value for the 
Markov Chain Monte Carlo (MCMC) analysis described 
below. Note that by fixing Tp, the MCMC code only 
needs to perform a full Keplerian orbit calculation once, 
rather than at every step in the chain, thus dramatically 
speeding up the computation. 

To estimate the uncertainty in the 10 remaining model 
parameters, we used an MCMC algorith m following 
the ba sic recipe described, for example, in I Winn et al.l 
(|2007f ). with whom we share the terminology used be- 
low. We produced 10 chains, each of 10^ points, starting 
from independent initial parameter values. The "jump 
function" was tuned so that ~ 25% of jumps were exe- 
cuted, and the first 10% of each chain was discarded to 
avoid the effect of the initial condition. We combined 
the chains, and in Table [T] we report the median value of 
each parameter, and assign uncertainties by taking the 
standard deviation of that parameter (except i, for which 
we report the 68.3% confidence interval). In Figure[21 we 
show the results of this MCMC analysis. One of the ma- 
jor benefits of MCMC is that probability distributions 
for derived quantities, such as the secondary eclipse im- 
pact parameter (in the bottom right panel of Figure [2|) , 
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Fig. 2. — Posterior probability distributions for Ri,, the planet- 
star radius ratio, Rp, and the secondary eclipse impact parameter, 
based on lO'^ MCMC samples. The secondary eclipse impact pa- 
rameter is the minimum projected separation between planet and 
star center (near time of superior conjunction), in units of the i?^. 
The dark shaded region shows where full secondary eclipses occur, 
while the light shaded region shows where grazing eclipses occur. 

are produced directly. Here, the impact parameter was 
taken to be the minimum projected separation between 
planet and star center in the units of the stellar radius. 

5. DISCUSSION 
5.1. Planetary mass, radius and density 
The planet's mass of 3.13 ± 0.11 Mjup ([Fischer et al.l 



l2007f ) is the third most massive among the known 
transiting planets, surpassed onl y by HD 147506b (8.6 
Mjup; also known as HAT-P-2b; I .Bakos et al.''2007^ and 
X0-3b (13.2 Mjup; iJohns-Krul l et al. 2007). We find 
that the mean density of HD 17156 is p — 3.8ti i g 
cm^"^, substantially higher than other transiting extra- 
solar planets, again with the excepti on of the remar kable 
planet HD 147506b (p=11.9 g cm-'^ HLoeillet et al.l2 007'). 
HD 17156b begins to bridge the gap between HD 147506b 
and the other transiting planets. 

The presence of a large rocky or metallic core has been 
suggested to account for "hot Jupite rs" with high densi- 
ties. F or HD 17156b, the models of iBodenheimer et al.l 
(|2003( ) predict a radius of 1.11 i?jup given the mass 
of HD 17156b, which is ~ 1 — a larger than the ra- 
dius derived from the light curve fitting, or cquivalently, 
the measured density is higher than the models would 
predict. However, for such a massive planet, this dis- 
crepancy is difficult to resolve by invoking the presence 
of a solid core, since the dependence of planetary ra- 
dius on the presence of a core is predic ted to be very 
weak in this rnass regime (e.g. Bodcnhei mer et al.ll2003l : 
Burrows et al.ll2 007: Fo rtnev et al.ll2007af l. For example. 



Bodenheimer et al., (2003) predict that the addition of a 
30 Mq core reduces the radius by only ^ 0.01 i?jup for a 
3 Mjup planet. This scenario would therefore require 
an extremely massive core to account for the measured 
radius of HD 17156b. 

Although a large rocky or metallic core can account 
for planets with high densities, the low density (i.e. 
large radius) planets such as TrES-4 ((Ma ndushev et al] 
[200% and WASP-1 ([Collier Cameron et al. 2007) seem 
to defy explanation. One possible mechanism for keep- 
ing planetary radii from shrinking during planetary evo- 
lution is to pump energy into the planet via tidal in- 
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Fig. 3. — Observed minus calculated times of transit for the four 
transit events discussed in this work. The dashed line shows the 
prediction assuming the values of Tc and P from Table[T] Since the 
value of P was determined from a simultaneous fit to the transit 
times and radial velocity data, it does not represent the best fit to 
the transit times alone. 



teractions (e.g. iMardlin^ I2007D . Such orbital energy 
transfer mechanisms depend strongly on the orbital ec- 
centricity, and HD 17156b represents a good test-case 
for these theories, having a significantly eccentric orbit 
(e — 0.67). However, the radius of HD 17156b is not 
in any way exceptional when compared to the existing 
extrasolar planets with close to circular orbits. This 
leads one to speculate that tidal energy transfer via or- 
bital dynamics may not play a substantial role in plane- 
tary radius evolution, or that in massive planets such as 
HD 17156b a mechanism is acting to allow the planets 
to more rapidly evolve to smaller radii. We note, how- 
ever , that the recently-anno unced massive planet XO- 
3b (| Johns-Krull et all 1200/1 also has an eccentric or- 
bit, but in contrast to HD 17156b, its radius is very 
large (1.92i?jup). Thus the challenges posed by newly- 
discovered exoplanets to theoretical models of their phys- 
ical structure seem to continue unabated. 

5.2. Search for transit timing variations 

The availability of multiple measured times-of-transit 
allows the presence of additional bodies in the HD 17156 
system to be detected or constrained, by searching for 
the influence of thei r gravitational perturbations on the 
orbit of HD 17156b (iHoIinan fc Murravl [2005t I Agol et all 
[200l ISteffen fc Ago]|l2005ll 

In Figure [31 we plot the observed minus calculated 
times of center of transit for the 4 events reported in the 
literature, including the event described in this paper. 
We find that these results do not yet reveal any evidence 
for timing variations that would indicate the presence of 
a third body in the system. Due to the long orbital pe- 
riod, HD 17156 may be more amenable to a search for 
such variations. However, the large orbital eccentricity 
implies a large region for which dynamical stability con- 
siderations would exclude the presence of such planets. 

5.3. Prospects for infra-red observations 

Our photometric analysis directly addresses the prob- 
ability that HD 17156b undergoes secondary eclipse, 
which is only possible for certain combinations of Ri, and 
i. We find a 9.2% chance for secondary eclipses to oc- 
cur, but only a 6.9 % chance for non-grazing eclipses. 
If, indeed, HD 17156b undergoes eclipses, i becomes 



locked rather tightly, thus greatly diminishing the al- 
lowed volume of parameter space for the stellar and plan- 
etary properties. Under this constraint, we find that 

the new best fitting values for {i,Ri,,Rp) are (88?7 qO^, 

1.34 ± 0.03 i?o, 0.89 ± 0.04 i?jup). These considerations 
serve to motivate a search for the secondary eclipses. 

HD 17156b's large eccentricity leads to a 25-fold vari- 
ation in received stellar flux over the course of the 
21.2169 d orbital period. This dramatic variation in illu- 
mination should produce complex weather at the planet's 
photosphere. In addition, the planet is subject to strong 
tidal forces during its periastron passage. The planet is 
therefore expected to be in pseudo-synchronous rotation, 
in which the planetary spin will be roughly synchronous 
with the orbit during the interval surrounding close ap- 
proach. The theory of Hut (1981) predicts a pseudo- 
synchronous spin period, Pgpin given by 



p . — 



1 f e2 



45 4 
8 ^ 



16" 



I + 3e2 



2)3/: 



r-Porbit = 91.3hr (4) 



Knowledge of Pspin, -^p, R-p, and the time-dependent 
pattern of received stellar flux make it possi- 
ble to compu te global climate mode ls for the 
plane t (e.g. ISho wman fc Guillot' '20 021: ICho et all 



Coooer & Showman 


2005, 


2007: iFortnev et al.l 


2007b: 



2003!; iBurkert et al. 2005 
2006; Langton & Laughlin 
Dobbs-Dixon fc Lin. ,2007, ). These models, in turn, 
can be used to obtain predictions of the planet's 
infrared photometric light curve at various wave- 
lengths. Model light curves can then be com- 
pared with obs ervations using the Spitzer Space Tele- 
scope (see, e.g.. ICharboiineau et"al]|2005t iDeming et al.l 



2005t iHarrington et alj l2006t iKnutson et al.l l2007allbl : 
Cowan et al.ll2007t ). 

We adopt the climate model of iLangton fc LaughlinI 
(|2007D and apply it to HD 17156b, using the known 
orbital and physical properties of the planet. The 
model uses a compressible 2D hydrodynamical solver 
(IM ams fc Swartztraubeij 119991 ) with a one-layer, two- 
frequency radiative transfer scheme. It is important to 
emphasize that our hydrodynamical model (like those of 
other workers in the field) is undergoing rapid develop- 
ment. Indeed, a primary goal for obtaining physical ob- 
servations of planets under time- varying irradiation con- 
ditions is to provide physical guidance to the numerical 
codes. 

Our model predictions are shown in Figure [H which, 
for reference, includes a model of the planet, the orbit, 
and the star plotted to scale (top left). The hydrody- 
namical model assumes a solar-composition H-He atmo- 
sphere for the planet, and is run for several orbits to 
achieve equilibrium, at which point we track the pre- 
dicted 8/um light curve for a full orbital period starting 
at apastron, by assuming a black body spectrum for each 
surface element in the model, and integrating the result- 
ing 8/im fiux over the visible hemisphere of the planet. 
The right column of Figure [4] shows a series of five global 
temperature plots, each corresponding to the thermal ap- 
pearance of the planet from the Earth. The figures are 
equally spaced in time by one-quarter of an orbit (127.3 
hours). 

The predicted 8/im flux from the planet is shown in 
the lower-left corner of the figure, and features a rapid 
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Fig. 4. — Upper left hand panel: The orbital geometry of 
HD 17156b. The small dots show the position of the planet at 
2.4 hour intervals throughout the 21.217 day orbit, and show the 
planet and star drawn to scale (ij* = 1.2 Rq). The locations of the 
planet at successive time intervals equal to one quarter of an orbit 
(1-5) are marked. The small orange spheres indicate the 91.3 hour 
predicted spin frequency of the planet, which is expected to be 
in pseudo synchronous rotation. Right hand panels: Temperature 
maps of the planet at locations 1-5, as viewed from the Earth. 
The temperature scale is mapped to a black body color map that 
approximates the planet's visual appearance in its own intrinsic 
radiation. Lower left panel: Predicted 8 fim light curve for the 
planet during the course of one planetary orbit. 

rise from a baseline of F-p/Fi, ~ 1.7 x 10^"* to F-p/Fi, ^ 
5 X 10^** during the ~ 30 hour interval following peri- 
astron. The magnitude of this rise is comparable with 



(but slightly less than) that expected for HAT-P-2b and 
HD 80606b (Langton & Laughlin 2007). 



6. CONCLUDING REMARKS 

HD 17156 thus represents a highly interesting system 
for follow-up studies in the mid infra-red (e.g. using 
the Spitzer Space Telescope) , even if it does not undergo 
secondary eclipses. Ground-based observations should 
also be pursued to continue the search for transit tim- 
ing variations, and hence to constrain the presence of 
additional bodies in the system, which, given the high 
orbital eccentricity of the planet HD 17156b, may show 
interesting dynamical properties and place constraints on 
planetary formation and evolution scenarios. This dis- 
covery represents an important success for the Transit- 
search. org project, and highlights the potential rewards 
of collaboration between distributed networks of amateur 
astronomers and the professional community. 
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